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We re-examine the kinematic variable mT2 and its relatives in the light of recent work by Cheng 
and Han. Their proof that mT2 admits an equivalent, but implicit, definition as the 'boundary of the 
region of parent and daughter masses that is kinematically consistent with the event hypothesis' is 
far-reaching in its consequences. We generalize their result both to simpler cases (ttit, the transverse 
mass) and to more complex cases (mTGen). We further note that it is possible to re-cast many 
existing and unpleasant proofs (e.g. those relating to the existence or properties of "kink" and 
"crease" structures in mT2) into almost trivial forms by using the alternative definition. Not only 
does this allow us to gain better understanding of those existing results, but it also allows us to 
write down new (and more or less explicit) definitions of (a) the variable that naturally generalizes 
mT2 to the case in which the parent or daughter particles are not identical, and (b) the inverses of 
ruT and mT2 - which may be useful if daughter masses are known and bounds on parent masses are 
required. We note the implications that these results may have for future matrix-element likelihood 
techniques. 



INTRODUCTION 

If a dark matter candidate is produced at the LHC, one 
of our experimental priorities will be to measure its mass. 
This is a non-trivial exercise, since dark matter is, by its 
very nature, invisible in the detectors. Thus, kinematic 
information is lost in each event. What is more, the 
difficulties are compounded by the fact that dark matter 
is invariably pair-produced. Despite receiving a lot of 
attention in the recent literature, we still only have one 
method for measuring masses^ in the case of pair decays 
which are identical, with each containing one invisible 
particle in the final state, ^ but are otherwise arbitrary 
[U |31 m |S].^ The method is based on the tot2 variable, 
introduced in |43l |44] and defined by 

171x2 = niin max(mT, m^). (1) 

Here mx and are the transverse mass variables for 
the individual decays, introduced originally for measur- 
ing the mass of the W^-boson [45 , 46 , and defined explic- 
itly below. The 'max' tells us to take the larger of these 
two variables. In the minimization, one is instructed to 
consider all possible partitions of the measured missing 



^ Once the dynamics are known or postulated, one could hope to 
measure masses directly from the matrix element. We comment 
on this later. 

^ For decays with more than one invisible particle, see 
^ In special cases, involving cascade decays and extra kinematic 
constraints, other methods are available '6', "T"!", '9', 'T0"TT"T2' 
[I3l[l4, 15, 16 17, 18 19 20, 21 22 23, 24 25 26, 27 28 29, 30 
[3Tl|32Ji33i^i35ii36ji37j^i39i|4Qii41j. An alternative method 
for the general case making use of the distribution of initial state 
radiation has also recently been proposed |42| . 



transverse momentum in the event between the two invis- 
ible particles, and to minimize with respect to partitions. 

This ad hoc definition is designed to cope with the 
fact that only the sum of the transverse momenta of the 
two invisible particles can be inferred from the missing 
transverse momentum observed in a collision, whilst in- 
heriting one desirable property from the usual transverse 
mass: it is bounded above by the mass of the parent 
particle. Unfortunately, mT2 is still not an observable 
when the masses of the invisible daughters are unknown 
and non-negligible, because the usual transverse mass is 
a function of the invisible daughter mass. In the case 
of neutrinos, whose mass can be neglected, this is not a 
problem. But it certainly is a problem for dark matter 
candidates. 

To cope with this obstacle, a futher ad hoc step was 
taken: consider as a function of the unknown invis- 
ible daughter mass m^: TOT2(wi). Now mT2{fni) is an 
observable, albeit an observable function. That is to say, 
each detector event returns a function. Unfortunately, 
in taking this step, the boundedness property of mT2 is 
lost: it is not true that mT2{mi) is bounded above by 
the mass of the parent, for abitrary values of mi. 

What then are the properties of mT2(f"'i)? On an 
event-by-event basis, mT2{'nii) is simply a smooth func- 
tion of nii . But if one plots the envelope of curves com- 
ing from many events, one discovers that the maximal 
curve features a kink [21 [31 IH [S]. That is to say, it is 
continuous, but not differentiable, exactly at the point 
•nii — rhi, where rhi is the true mass of the invisible 
daughter. Moreover, since we already know that the 
maximal value of 771^2 (w^) is the (true) parent mass, toq, 
we see that the kink has co-ordinates {rfii, mo); by identi- 
fying the location of the kink in an experiment, one may 
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measure the masses of both the parent and the invisible 
daughter. 

On reflection, this result, though pleasing, is somewhat 
mysterious. One started from an ad hoc definition of a 
transverse mass for identical pair decays, guaranteeing 
only the desirable property that 771^2(^1) ^ wq- More- 
over, generalization to rrii ^ rhi, in which the one de- 
sirable property contained in the definition is lost, gives 
rise to a strange kink behaviour. What does all of this 
mean? 

Recently, Cheng and Han, gave an elegant interpreta- 
tion of the function mT2{'mi) llZj-'* They showed that, 
for a given event, it defines the boundary of the region 
in the (771^,7710) plane for which the various kinematic 
constraints, namely conservation of four-momentum and 
the mass shell constraints, admit a solution. By 'admit 
a solution', one means that there exist real values of the 
unknown momenta, and real, non-negative values of the 
unknown energies, solving the constraints. The existence 
of such a solution for a given value of (777^, 777.0) means that 
one cannot rule out the possibility that the true mass val- 
ues, (777^,7710), are given by (777^,777,0), on the basis of the 
information obtained from that event. 

The proof is very simple, though we refrain from re- 
peating it here (we shall, in any case, give a proof for 
a more general case of non-identical pair decays in what 
follows) . The beauty of the result is that it shows that 
the original ad hoc definition of 7777^2 ("^i)i and its ad hoc 
extension to 7777^2 (ttt.^), fortuitously give rise to a natural 
function, namely the function that defines the boundary 
of the allowed region in mass space, on an event-by-event 
basis. If one considers multiple events, the allowed re- 
gion is restricted to the intersection of the allowed regions 
coming from each event. If one considers arbitrarily many 
events, one ends up with an extremal allowed region, and 
a corresponding extremal boundary. In the case of iden- 
tical pair decays this extremal boundary is precisely that 
containing the kink identified in [21 121 IH |S] . 

In what follows, we would like to show that this re- 
interpretation of 777^2(777^) in terms of a boundary is both 
general and powerful. Firstly, we remark that it also ap- 
plies to the simpler case of the usual transverse mass 
variable for single-particle decays. Secondly, we show 
how it can be used to give an almost trivial derivation 
of the form of the kink curve for identical pair decays. 
Thirdly, we show how it can be generalized to pair de- 
cays in which either the parents, or the invisible daugh- 
ters, or both, are not identical, and have different masses. 
This generalization may be practically useful, for exam- 
ple, in the case of squark-gluino production in the con- 
text of supersymmetry, or in theories in which more than 
one particle is stable on the length scale of a detector. 



* A similar interpretation was given without proof in |30| . 



In these cases, the number of masses that are, a pri- 
ori, unknown is increased. Consequently, the mass space 
is higher-dimensional, and so is the boundary of the al- 
lowed region that follows from applying the kinematic 
constraints to an event. Nevertheless, the form of the 
extremal boundary is easily obtained. For the case of 
distinct parents (such as a squark and a gluino) decay- 
ing to a common LSP, the extremal boundary forms a 
surface in the three-dimensional space parametrized by 
the three unknown masses. The surface is, as we shall 
see, creased, with various kink structures visible in two- 
dimensional projections. Fourthly, we explain how ob- 
servables introduced previously to cope with combina- 
torics and upstream or initial state radiation can also be 
understood as generating the kinematic boundary. 

We stress that our arguments are purely theoreti- 
cal, and take no account of what might realistically be 
achieved in kinematic measurements at the LHC. Nev- 
ertheless, we feel that, in order to properly understand 
and use variables like 771t'2, it is important to know both 
how the variables behave under ideal conditions, as well 
as how this is modified in real situations. Although we 
shall not address the latter aspect here, we hope to make 
a useful contribution to the former. 

As an example of how reality deviates from the ideal, it 
is clear that a real LHC data sample cannot saturate the 
true extremal boundary. Indeed, the extremal boundary 
corresponds to events in which the parent particles are 
infinitely boosted with respect to the laboratory frame [5] 
by radiation upstream or in the initial state. Neverthe- 
less, even a subset of events, such as those contained in 
a finite LHC data sample, defines a corresponding kine- 
matic boundary. The arguments we give show that re- 
construction of that boundary is the best that one can 
hope to achieve in the absence of additional kinematic or 
dynamic information, whether inferred or assumed. As 
another example, it is not yet clear how the boundary hy- 
persurface we describe might best be reconstructed from 
LHC data. What one would like to do is to generate one- 
dimensional distributions of observables, which can then 
be fitted by Monte Carlo simulations. Towards the end 
of the paper we make a partial effort to address this, by 
discussing how such observables, related to the boundary 
hypersurface, may be derived. 

Our notation is as follows. For a single particle de- 
cay, we consider a parent particle of mass 7770 decay- 
ing into an invisible daughter particle of mass 777^ and 
a system of visible daughter particles of invariant mass 
niy. We write the four-momenta of particle by Po = 
(i?o,Poj 9o), where po is transverse to the beam direction 
and go is parallel to it. We denote the transverse energy- 
momentum by ao = (eoj Po)i where the transverse energy 
is defined by eg = Po + "^o- ^'^^ P^i^' decays, we use un- 
primed quantities for one decay and primed quantities for 
the other. We shall often need to distinguish between hy- 
pothesized values of the unknown masses of the parent 
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and invisible daughter, and the true values; we denote 
the latter with tildes. 

Finally, to avoid confusion, we remark that we shall al- 
ways illustrate our arguments with the special case where 
the visible daughter system contains a single, massless 
particle. Although theoretically the simplest case, this 
is probably the least favourable example from an exper- 
imentalist's viewpoint, since a kink is generated in this 
case only by events in which there is significant upstream 
transverse momentum. Nevertheless, our general argu- 
ments apply to arbitrary visible systems, including those 
which appear to be experimentally more favourable. 



SINGLE PARTICLE DECAYS AND THE 
TRANSVERSE MASS 

Let us first prove that for an event consisting of a single 
particle decay, the locus of the curve mp — mximi) is 
equivalent to the boundary of the region in (m^jTOo) for 
which the kinematic constraints^ 
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Po = {Pi+Pvr = m 



2 _ ^2 
0' 



(2) 
(3) 
(4) 



admit a solution, in the sense defined above (with real 
momenta and real, non- negative energies). In the above, 
and ^ are measured, whereas Pi are four unknowns. 
Here, mx is defined by 



rriy = (a„ -|- aif 



2(e.„ei - p„ • Pj). (5) 



The proof is in two parts. First, we establish that 
any (m^, mo) for which (2|4) have a solution is such that 
Too > rnximi). Second, we establish that (2]|4| have a 
solution for toq = mT{mi). 

For the first part, we have 



ml = {pv + Pif ^ pI + pj + "^-Pv ■ Pi (6) 
2{EyE, - p„ • Pi - q^qi). (7) 



2 I 2 

TO,, + TO, 



But since E^Ei — q^qi > e^Ci (with equality at E^qi — 
Eiqv), we have that {a^ + onY < [Pv + PiY, or, in other 
words, mT^rrii) < too. 

For the second part, we need to show that the equa- 
tions 
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Pt =m,, 

2 _ ™2 f^2\ 



{Pi+Pv) =TOj,(toJ, 



(8) 
(9) 
(10) 



admit a solution. Equations (10 1 fix pi, and equation 
([s]) fixes Ei in terms of qi. Equation ([9| is satisfied by 
requiring {pi+PvY = {av + otif' , which, as we just learnt, 
requires = j^. Now for a given event, takes a 
value in [—1,1], and as the remaining unknown qi varies 
in K, ^ takes all values in [—1,1]. Thus, equations (8|l0l 



admit a solution 

So for a given event, too — mT{mi) defines the bound- 
ary of the allowed region in (mi, mo). Given multiple 
events, the allowed region shrinks to the intersection of 
the allowed regions for each event. In the limit of ar- 
bitrarily many events, we obtain the extremal boundary 
given by the kink curve. Its explicit form was derived in 
[31 H]. Here we simply quote the result for the special 
case where the visible system consists of a single, mass- 
less particle. The locus of the extremal boundary is given 

by 

for TO 
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TO,; > TO, 



(11) 

(12) 



These are simply straight lines in the space of mass- 
squareds. 



IDENTICAL PAIR DECAYS AND m.T2 

The analogue for identical pair decays and TOr2 of the 
argument just given has already been given in [47], and 
in any case follows as a corollary from our analysis of 
non-identical pair decays to be given below. Here, we 
point out that it can be immediately used to derive the 
form of the extremal boundary or kink curve, from the 
form of the maximal curve for a single decay (given by 
( TI|T2 | for the special case of massless visible particles). 

Indeed, for a single event, the m^i locus is given by 
the boundary of the region in (mi,mg) for which the 
equations 
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(13) 
(14) 
(15) 
(16) 
(17) 

have a solution, with to^ = to^ and TOq = too. The ex- 
tremal boundary is given by the intersection of the al- 
lowed regions for all possible event configurations. But 
in considering all events, we permit any value of j4, such 
that the last equation (17 1 can always be satisfied and 
becomes trivial.^ But then the remaining equations de- 
couple into those for the individual decaying systems. For 



Cheng and Han |47| call these the minimal kinematic constraints, 
because one may have supplementary constraints in theories with 

cascade decays. ® This does not hold if wc only consider a subset of events, for 
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FIG. 1: Representation of the bounding planes (visible faces) 
and the extremal allowed region (solid) for the case described 
in the text with rhi = fh'i, rrii = rn'^, and m„ = = 0. The 
vertex representing the true values of the masses is indicated 
with a red ball. The origin of the axes is at the point (ttiq = 
jfiQ — rhi,m'o = ifiQ — , = 0). 



these systems, the extremal curve has ah'eady been com- 
puted in Since we have identical decays, the extremal 
curves for the two individual systems are the same, and 
this same curve is the extremal curve for the pair de- 
cay. By this argument, which we call the 'decoupling 
argument', we thus find a simple proof of the results pre- 
viously obtained in [H |5] . 



NON-IDENTICAL PAIR DECAYS. 

The most general case is the one in which neither the 
parents nor the invisible daughters have common mass. 
The kinematic constraints are then given by ( TspT I, but 
now with ^ rrii and mg ^ mo- It is simple enough to 
derive the form of the extremal boundary of the allowed 
region in such a case. Just as for identical pair decays, 
we may invoke the decoupling argument. The extremal 
locus is then given simply by the individual extremal loci 
for the individual decays. Thus, in the four-dimensional 
space (mi, toJ, mo, mg), the extremal locus (which is a 
two-dimensional surface) simply factorizes into the prod- 
uct of the two one-dimensional curves in (m^jmg) and 
(m^,mg). Things become more interesting if we posit 
that either the parents or the daughters have a com- 
mon mass. For example, if there exists a pair-produced, 
unique, stable dark-matter candidate, then the daughter 
particles in a pair decay will have common mass. Let us. 



example those in which initial state or upstream momentum is 
forbidden, such that p + + p'^^ =0. 



for the sake of argument, consider this case in what fol- 
lows. The allowed regions now occupy the 3-dimensional 
space parametrized by (mo,mQ,mi = m[). What is the 
form of the extremal boundary surface? By the decou- 
pling argument, the extremal allowed region is given by 
the intersection of the extremal allowed regions for the 
individual decays. In the special case of massless visible 
particles, for example, the individual extremal allowed 
regions are bounded by 



ml 



™2 



m. 



ml mo < mo, 



m. 

-r^, mg > mo, 
mf 



for the unprimed system, and by 



mn 



m?, 



niQ < rfiQ 



(18) 
(19) 

(20) 
(21) 



for the primed system. These both describe a surface in 
the space with co-ordinates (mg, m'^^mf = m'^). The ex- 
tremal boundary for the pair decay is then given at each 
point (mQ,mg^), by the surface that gives the smaller 
value of mf. Let us assume that mg > mg, without loss 
of generality. The extremal boundary surface is sketched 
in Figure [T] The intersections of the four planes give rise 
to four creases in the extremal boundary surface, gener- 
alizing the kinks observed in extremal curves for identi- 
cal pair decays. Furthermore, one can see that various 
types of kink behaviour may arise by taking various two- 
dimensional slices through the three-dimensional space 
of masses. 

Consider, for example, the extremal curves of mg vs. 
mi obtained at fixed mg. For mg > mg there are two 
kinks - the 'usual one' at (mg = mo,mi = rhi) and a 
second on the upper part of the diagonal crease in Figure 
[l] for m§ — mf < mg < mg there is a single kink on the 
lower part of the diagonal crease; for mg < mg — mf the 
allowed region is null. The locus of the diagonal crease 
in (mQ,mg^) is given by 



™'2 



mn - mn 



"0 



ml, 



~ 2 ~ 2 ^ 2^-2 

mr, — m, < mn < m 



0' 



mo > m. 



^0- 



(22) 
(23) 



A particularly striking kink is seen if one fixes mf and 
considers mg^ as a function of mg. There is then a right- 
angled kink coming from the diagonal crease. 

In the above, we have chosen to derive the extremal 
boundary directly, rather than relying on an explicit 
functional definition (like the original definition of mT2)- 
Nevertheless, an explicit functional definition is easily 
guessed. If we just have non- identical daughters, but 
common mass parents, we can simply use the original 
mT2 definition. But if the parents have distinct mass. 
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the original tot2 is no good, because its definition ex- 
ploits the equality of the maximal values of rriT and mip. 
To maintain this equality, we propose a generalized defi- 
nition of tot2 for non-identical parent decays as 



m|2(rnj,m-, Too/too) 



mm max 



mo 



mo 



m'n 



(24) 



Note that both ^m^^fhi) and ^m'^{rh^) are bounded 
above by momQ. 

It remains to show that this definition reproduces the 
boundary of the allowed region, event by event. As 
for single particle decays, there are two parts to the 
proof, which generalizes immediately from that given 
in [47j. First, we establish that any (rrii, m^, mo, TOq) 
for which ( 13][T7 l have a solution is such that mom'n > 



m|n2(wi, m^, mo/mo). This follows immediately from 
m|n(mi) = (a„ -|- a^)^ < (p^ + Pi)'^ = mg, from the cor- 
responding inequality in the primed system, and from 
m|,2 < max(^m|, ^m^) < momj,. 



Second, we establish that ( 13p7 l 



tion for 



m, 







^m|2(»7ij, m-, mo/mo) 



have a 
and 



solu- 



^m|,2(mi, m^, mp/mo). There are three possibilities to 
consider, arising from the three different ways in which 
values of mT2 niay arise [H]: (i) the balanced case, with 

(ii) the unbalanced case with 



'■T2 



^T2 



> ^mjn, and (iii) the unbalanced case 
i'^. In case (i), the solution of 



> 



with m, 

( 13][T7 ) is given by the pi assigned by the minimization 



in the definition of mT2 aud the qi such that — 



1^ 

_Bi ~ _E„ ' 

and similarly for the primed quantities. In case (ii), the 
solution has — -g^; to find a suitable q^, we note 
that, if we chose qi such that j^- — j^, we would ob- 
tain {pi + PvY — ™T < ^"^T2' whereas if we chose 
qi oo, we would find {pi + PvY ~^ Since {pi +PvY 
is a continuous function of qi on K, there must, by the 
intermediate value theorem, exist values of qi such that 
{Pi + Pv)^ = ^'m'T2^ ^ required. A similar argument 
applies to case (iii). 

Thus we have proven that the explicit definition of 



mT2 for non- identical decays in ( 24 1 reproduces the kine- 
matic boundary surface. Interestingly, it would ap- 
pear that, in the case of distinct daughter and par- 
ent masses, the boundary for each event is given by a 
three-dimensional hypersurface in the four-dimensional 
space of (mo, mjj, mi, m^), whose locus is momj, = 

m|i2('7ti, m^, ^)- By contrast, the extremal houndaiy in 
this case, whose form we derived at the beginning of this 
Section, is a two-dimensional surface in four dimensions. 



COMBINATORICS AND rriTGc 



In the real world of experiment, one must also face the 
fact that in pair decays there will be combinatoric ambi- 
guities. For example, in identical pair decays, there are 
(at least) two copies of each visible particle in the final 
state, and one does not know which decay to assign them 
to. With 2n visible particles, this results in a 2(»-i)-fold 
ambiguity. As a result, the visible momenta in equations 
( TsfrZ l are themselves ambiguous. In the presence of 
such an ambiguity, the allowed kinematic region is given 
by the union of the allowed regions obtained by consid- 
ering all possible branch assignments. Equivalently, the 
boundary of the allowed region is obtained by taking the 
minimal my2 curve with respect to the different com- 
binatoric assignments. This corresponds to the niTGen 
variable defined in 

Similar combinatoric ambiguities can arise in the pres- 
ence of upstream or initial state radiation. Again, to 
find the allowed kinematic region, one simply takes the 
union of regions obtained by considering all possible as- 
signments, which is equivalent to the mT2-based prescrip- 
tion given in ,49] . 



DERIVED OBSERVABLES AND 
ONE-DIMENSIONAL DISTRIBUTIONS 



The kinematic boundary hypersurface for a decay 
topology gives a complete picture of the kinematic con- 
straints coming from an event. In principle, it can be 
generated from experimental data, but in practice, it is 
not clear how this will be achieved. More likely is that ex- 
perimentalists would prefer to generate one-dimensional 
distributions of specific observables, which they can then 
compare with the results of numerical simulations. 

In order to do so, one would like to understand how to 
'translate' the kinematic boundary plot into observables. 
This can certainly be achieved if one has determined all 
but one of the unknown masses. One simply needs to iso- 
late the pertinent observable. For example, in an iden- 
tical pair decay, if one knows the mass of the invisible 
daughter, the mT2(mi) distribution can be used to mea- 
sure the mass of the parent. Conversely, one could imag- 
ine that one knew the mass of the parent, and wished to 
extract the mass of the invisible daughter. Clearly what 
one needs in this case is an observable derived from the 
inverse function of 771x2(1^1) ■ 

To derive an explicit expression for the inverse function 
of mT2i consider first a single decay and my. We claim 
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FIG. 2: Distribution of the usual mT2 for mo = 400, 
m-o = 600, rrii = 100, i.e. for a case in which the parent 
particles have different masses. The distribution has a kine- 
matic endpoint at the mass of the heavier particle (indicated 
by the red arrow) but it will be noted that this endpoint is 
not very useful, as the density of states near the end point is 
low. 



that the inverse of rax for an event is given by'' 
(m^)~^(mo) = (ao - a^Y , TOq > m.^. 



(25) 



(The condition toq > is added to ensure that mT{rai) 
is surjective, such that exists.) The reader may 

easily check expHcitly that 771^^(771^(^1)) — ™i- 
For pair decays, the appropriate definition is 



''T2 



2(7770) = max min(777j,^, 777^ ^). 



(26) 



Note that, rather than minimizing over momentum as- 
signments, we now maximize over them, and, rather than 
taking the larger of the two observables, we now take the 
smaller of the two. To demonstrate explicitly that this 
is the inverse, it suffices to show that both 771^2 and its 
inverse give rise to the same kinematic boundary for an 
event. 

We note that the inverse of 771t2 can be applied directly 
to measure the mass of a common invisible daughter, 
even if the parent masses are not identical, provided the 
parent masses are themselves known. 

Finally, let us ask what are the appropriate observ- 
ables for non-identical decays? If one knew, say, the two 



^ We note that the inverse of rriT may also be useful on its own 
for hadron collider mass measurements. For example, one could 
measure the charged Higgs mass in decays t — > H'^b, followed 
by the decay — ► tv with an invisible daughter. In fact, the 
variable defined in [50) is precisely the inverse transverse mass of 
the tbH~^ system, with parent t and invisible daughter . 



daughter masses as well as the ratio of the two parent 
masses, one could simply invoke the observable defined 
in ( 24 ) to determine the product of the parent masses. If, 



alternatively, one knew the two daughter masses as well 
as the lighter parent mass, one could use the usual 7Mt2 
definition. Indeed, it is easy enough to see that 777^2 is 
bounded above by the mass of the heavier parent. Unfor- 
tunately, the 777T2 distribution has a very poor endpoint 
behaviour, as we illustrate in Figure [2] Experimentally, 
such fine edges are likely to be difficult to catch. 



Since the distribution of the ratio variable (24 1 has 
much better end-point behaviour - it has thick edges, 
especially when the correct value of the ratio is used - the 
pragmatic choice may be to hypothesise different values 
for the ratio, and to com par e the experimental data to 



template distributions of ( 24 1 for each value of ^ . 



MASS MEASUREMENT VIA THE MATRIX 
ELEMENT 

If one knew (or guessed) explicitly the Lagrangian, one 
could hope to measure the masses directly using the like- 
lihood and the matrix element. In practice this reduces 
to kinematic constraints plus constraints from parton dis- 
tribution functions (PDFs), plus a small dependence on 
dynamics (spins/couplings etc), which in any case one 
would presumably not want to trust to begin with. The 
arguments of Cheng and Han tell us that all of the infor- 
mation from kinematics is encoded in 7777^2, which cuts off 
the lower right region of the {mi,mo) plane. The PDFs 
limit the allowed mass of the parent, cutting off the upper 
part of the plane. So one can already see, roughly speak- 
ing, what the negative likelihood contours of a matrix 
element method will look like. They will describe a nar- 
row gully lying along the boundary of the extremal 771T2 
curves. Evidence is beginning to emerge which supports 
this conjecture |5T| . 



SUMMARY 

Cheng and Han's interpretation of the function 77772 
as the kinematic boundary between allowed and disal- 
lowed regions of mass space is a powerful one. It has 
allowed us to prove more elegantly the known results for 
the extremal boundaries for single-particle and identical 
pair decays. Its generalization has allowed us to prove 
new results for non-identical pair decays and for complex 
pair-decay topologies with indistinguishable particles in 
the final state. 

We have constructed three explicit examples of bound- 
ing functions that perform roles similar to 777^ and 77772 > 
but with differing assumptions. The first is a general- 
ization of 77772 that is appropriate when parents with 
different masses decay to equal-mass invisible daughters 
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- a case which will be of particular interest at the LHC. 
The other generalizations are the inverse functions m^^ 
and which require the parent particle mass as a pa- 
rameter, and which then provide the extremal bound on 
the invisible daughter (WIMP) mass. 

As well as providing a mass-determination method in 
their own right, such variables encode the kinematic part 
of the likelihood function. This means that insights 
gained from their construction can inform one's inter- 
pretation of mass determinations using the full matrix 
element - where such calculations are computationally 
tractable. This final test will show whether it is safe to 
neglect the effects of spin, determine the character of the 
creases, and get the desired results by using the bound- 
ary. 
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